Generalized fluctuation-dissipation relation and eff^ective temperature 

in off'-equilibrium colloids 
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The fluctuation- dissipation relation (FDR), a fundamental result of equihbrium statistical physics, 
ceases to be valid when a system is taken out of the equilibrium. A generalization of FDR has been 
theoretically proposed for out-of-equilibrium systems: the kinetic temperature entering FDR is 
substituted by a time-scale dependent effective temperature. We combine the measurements of 
the correlation function of the rotational dynamics of colloidal particles obtained via dynamic light 
scattering with those of the birefringence response to study the generalized FDR in an off-equilibrium 
clay suspension undergoing aging. We i) find that FDR is strongly violated in the early stage of the 
aging process and is gradually recovered as the aging time increases and, ii), we determine the aging 
time evolution of the effective temperature, giving support to the proposed generahzation scheme. 



The study of the dynamics in non-equilibrium systems 
is an intriguing and fascinating area of modern physics. 
There is an obvious interest in the off-equilibrium regime 
as condensed matter is often found far from ideal equilib- 
rium condition. In these cases where "things keep hap- 
pening on all time scales" Q some of the most funda- 
mental symmetries of equilibrium statistical mechanics 
are broken. A particular class of systems where the out- 
of-equilibrium status occurs naturally are the so called 
glass forming systems. These systems include, for ex- 
ample, disordered spin systems close to the spin-glass 
transition, supercooled molecular liquids, and jamming 
colloidal solutions. In general the dynamics of these sys- 
tems is extremely sensitive to change of external param- 
eters so that, when they are a little cooled or densified, 
the time required to re-establish equilibrium grows enor- 
mously and an off-equilibrium regime is entered. Here 
the average quantities becomes time dependent and the 
correlation and response functions depend on two times, 
in this regime the system is known to perform aging). 

Correlation functions and response functions are the 
basic quantities through which one probes the dynamics 
of a system in many-body statistical physics. These func- 
tions are closely related in equilibrium by the fluctuation- 
dissipation theorem (FDT). The FDT Q establishes a 
relationship between the correlation function CAsit) — 
{A{0)B{t)) of the spontaneous equilibrium fluctuations 
of the dynamic variables A and B and the response func- 
tion XAB{t) = {A{t))/h, describing the change in the 
average value of A due to an infinitesimal external field 
h that is coupled to the variable B in the perturbation 
Hamiltonian: 



XAB{t) = f3[CABiO) ~ CAsit)] . 



(1) 



Here (3 = l/kgT, where T is the kinetic temperature of 
the system and A and B are variables that have zero- 
mean in the unperturbed case. In this formulation of 
the FDT the field h introduces an energy contribution 
6H = —hB in the system Hamiltonian and it is switched 
on instantaneously aX t — Q and kept constant for t > 
as an Heaviside step-function. In the following we will 
refer to auto-correlation function [B = A) and drop the 
label A A in CAA{t) and XAA{t)- 

Out of equilibrium the system is non- stationary and 
time-translational invariance is lost. The correlation and 
the response become two-times quantities depending also 
on the aging time tu,'. C = C{tyj, tyj-\-t), x = X^w^tw^f) 
and the FDT ([T]) is not supposed to hold. The impor- 
tance of extending the theorem to the non-equilibrium 
regime has led to propose the generalized fluctuation- 
dissipation relation (GFDR) 3, 4, 5]. 

The generalization of the FDT proposed by Kurchan 
and Cugliandolo @, 0] in the early '90 can be expressed 
as follows: 



x{iw,tw + t) ^ [3 



dsX{t^,.s)dsC{t^,s) (2) 



^Electronic address: cmaggi@ruc.dk 



where for short times t (i.e. t/t^ ^ 1) the previous 
equation reduces to the FDT (Eq. [T]) and the system is 
said to be in a quasi- equilibrium state; for intermediate 
t (tjtw — 1) the violation function Xlt^^t) quantifies 
the deviation of the GFDR from the FDT; finally, when 
t/tyj ^ 1, the function X depends on i^, and s only 
through the function C. The violation function can also 
be interpreted 0, 0] in terms of an effective temperature 
Teff = T/X. The latter has complicated behavior for 
intermediate t values, but in the long-i regime it assumes 
a value that depends only on the waiting time t^ and 
one can think of it as "waiting-time dependent effective 
temperature" . 

The GFDR in out-of-equilibrium systems has been 
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studied theoretically and through computer simula- 
tions in spin-glasses [13, El and in models for glassy 
dynamics [12{. Later, off-equilibrium molecular dynam- 
ics simulations have given the possibilit y to perform the 
same analysis on structural glasses 13, Re- 
cently the generalized relation has also been investigated 
in more exotic systems, like in the simulations for a glassy 
protein [l7|, in self-assembling processes of viral capsids 
formation and of sticky disks crystallization [3] and in 
active systems composed of self-propelled particles [l^l ■ 

On the other hand few experiment have attempted to 
study the GFDR[23, HI 0, [13, S [H ■ Furthermore 
the results reported until now for structural glasses re- 
main controversial [H HI, [H, [H [H • The experiments 
in structural glasses are intrinsically difficult because it is 
necessary to simultaneously measure a correlation func- 
tion of a given variable and the associated response func- 
tion, and all this in a systems that is instantaneously 
brought out-of equilibrium. These difficulties are relaxed 
in the case of colloidal glasses (or gels, or jams), because 
the associated timescales are much longer and fall easily 
in the experimentally accessible window. 

In this letter we report an experimental investigation of 
the generalized FDR in a colloidal systems composed of a 
water suspension of clay platelets (Laponite), which is off 
equilibrium when it ages towards the final arrested state 
(gel or glass, depending on the clay concentration [13] )■ 
We study the reorientational dynamics of the asymmet- 
ric clay platelets, looking at the orientational correlation 
functions via depolarized dynamic light scattering and at 
the corresponding response function via the electric field 
induced birefringence. Measuring both C = C{tw,tw+t) 
and X = xi^w, tw +i) at different waiting time during the 
(days long) aging process we find that FDR is strongly vi- 
olated in the early stage of the aging process and is grad- 
ually recovered as the aging time increases. Moreover, 
from the parametric C— x plot (the so called "FDT plot") 
we determine the effective temperature and follow its evo- 
lution from the high values (Te///T « 5) pertaining to 
young systems towards the equilibrium [T^ff/T = 1) at- 
tained at long waiting time. Our findings confirm the 
generalization of the FDT to off-equilibrium systems pro- 
posed by Cugliandolo and Kurchan 15 year ago. 

The off- equilibrium sample. A solution, prepared by 
stirring the Laponite powder with water, evolves toward 
an arrested state on a timescale that span the ranges 
hours to months when simply kept at room temperature 
and pressure. Even low concentration aqueous solutions 
of this colloid, as the one used in our experiment (~ 1 % 
Laponite weight fraction) , show strong aging of its light- 
scattering correlation function [sO] ■ Due to its long (with 
respect to the experimental timescale and to the decorre- 
lation time) aging process of the systems, one can approx- 
imate the different measure of fluctuations and response 
as obtained in a sequence of steady out of equilibrium 
states. Furthermore, the anisotropic shape of the clay 



disc makes it possible to study its reorientation dynam- 
ics through the response and the correlation function. 
Laponite particles dissolved in water have the form of 
fiat cylinders with a diameter of 25 nm and an thickness 
of about 1 nm. Laponite colloidal particles are good scat- 
terers of visible light and this allow us to rapidly measure 
the autocorrelation function of the scattered field. 

The correlation function. In a dynamic light scatter- 
ing experiment one measures the correlation function of 
the optical field scattered by the sample. The scattered 
field can be directly related to the translational and ro- 
tational motion of the anisotropic colloids suspended in 
the solvent [13]. The colloid's rotations are related to the 
second rank tensor of the optical susceptibility. Specif- 
ically, in the VH (depolarized) scattering geometry, one 
measures the autocorrelation function of a variable that 
depends on the platelet 's orientation: 



A{t) = Y,P^{cos{e,{t))) 



(3) 



where P2{x) — [ix^ — l)/2 is the 2nd order Legendre 
polynomial, 9i is the angle formed by the symmetry axes 
of the i-th particle with the polarization vector of the 
incident field and the sum is extended over the parti- 
cles contained in the scattering volume [l^l ■ This results 
holds exactly only if the time scale of the rotational dy- 
namics is much faster than the translational one. This 
assumption was confirmed by comparing the VV (polar- 
ized) and VH (depolarized) photon correlation (PCS) at 
different waiting times and clay concentration [29|. The 
autocorrelation function of A was measured using the 
VH geometry via PCS. Several autocorrelation functions 
have been measured during the aging process of the sam- 
ple with a time resolution (1 /xs) dictated by the time- 
structure of the detector (photomultiplier) response. 

The response function. If one applies an external field 
that tend to align the particle the system -due to the 
anisotropic platelet 's polarizability- becomes birefringent 
31, 111- If the aligning field is a DC (or low frequency) 
electric field (Kerr effect) the degree of rotation of a lin- 
early polarized laser beam is proportional to the square 
of the electric field via a coefficient that is proportional 
to A (Eq. [3|). Therefore, the (time dependent) Kerr re- 
sponse to the switch-on of an electric field is proportional 
to the desired response function (i.e. the response conju- 
gated to the correlation function measured in depolarized 
PCS). For selected values of the waiting times, the time 
resolved response functions and the corresponding corre- 
lation functions, was measured during the aging process 
of the Laponite solution. The length of the electric pulses 
produced sets the dynamic window of our experiment to 
about 1 ms. 

Note that the relaxation time r of these functions 
is always much smaller that the typical waiting time 
(r ^ t^). This means that the time-resolved correla- 
tion and response are well-defined quantities although 



3 



PMT 



CORRELATOR/ 
COMPUTER 



COLLIMATOR- 



SCATTERED LIGHT- 



PS 

=> L3 



LASER 



P1 



WINDOW- 
KERR CELL- 
ELECTRODE- 



QWP P2 PHOTODIODE 

H I 



' L2 



ELECTRIC 
PULSE 



AMPLIFIED 

PULSES 

SOURCE 



DIGITAL 

OSCILLOSCOPE 



FIG. 1: Sketch of the experimental set-up (see Ref. [29| for more 
details). The laser's radiation (A = 532 nm) is polarized by the 
polarizer PI and focused by the lens LI at the centre of the cell 
containing the sample. The scattered light is collected by the lens 
L3 and selected by the polarizer P3 (orthogonal to PI). A photo- 
multiplier tube (PMT) detect the scatterd photons. When no elec- 
tric pulse is applied to the cell the output of the PMT is acquired 
by a computer equipped with a custom digital correlator, this mea- 
sures and stores the correlation function. The Kerr-cell containing 
the sample is provided by two electrodes connected to a source of 
amplified electric pulses. The forward-scatterd light (rotated by 
the electrically stimulated sample) is collected by the lens L2 and 
passes through a quarter-wave plate (QWP) and the polarizer P2 
(orthogonal to PI). The transmitted light is detected by a photo- 
diode connected to a digital oscilloscope. This is triggered to the 
source of electric pulse measuring and storing the Kerr-response 
function. 



the system is aging. In addition, if any FDT-violation 
can be detected, we expect to find that in on a timescale 
comparable to the relaxation time [t/r > 1). 

Examples of the measured quantities are shown in the 
left panel of Fig. [21 The correlation function and the 
response function are reported as functions of t for dif- 
ferent aging times t^. For short t the FDT holds while 
we can see a clear deviation from the FD relation for long 
t where x/f3 does not overlap with 1 — C {C{t) is normal- 
ized to C(0) = 1). When x{tw,tw+t)/ (3 is parametrically 
plotted against C (t^ ,1^^ +t) using t as parameter (FDT 
plot) the departure from the 1 — C line becomes evident 
(see the right panel of Fig. [2|). The deviation from the 
behavior expected from the FDT reduces its importance 
as grows, and the time where Tx and C detaches from 
each other moves to longer t (see also Fig.[31Ja), where the 
interested region of the FDT plot has been expanded). 
In order to quantify this deviation, we have performed 
a straight line fit to the longer time points in the FDT 
plot. The slopes (m) of these lines are a measure of the 
effective temperature: l/m = Teff/T. 

The tm dependence of T^fj is reported in Fig. [3Kb): 
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FIG. 2: (Left) Normalized time correlation function (o) and re- 
sponse function (□) measured at different aging times (from top to 
bottom: tm = 90, 125, 225 and 1200 min), the solid line represents 
l — C. Note that 1 — C deviates from x/P when the FDT is violated. 
The importance of this deviation reduces as increases. (Right) 
Response function vs correlation function measured at different ag- 
ing times (o), the black line represents the expectation of the FDT 
while the color lines represent the linear fits to the points in the 
off-equilibrium section of the FDT plot. It can be appreciated 
how these points corresponding to long times approach gradually 
the FDT as the aging time increases, for the longest waiting time 
the fitting line overlaps almost perfectly with the prediction of the 
FDT (see Fig. \3\ for a comparison of the FDT plots at different 
aging times). 



Teff decreases as increases. The linear fit to the long 
t region of the FDT plot also defines a characteristic 
value of the correlation C where the FDR breaks-down, 
the so called Edwards- Anderson value q; this quantity is 
reported as a function of in Fig. [Sjc). Finally, the 
quantity q, via C{to) — q, identifies a characteristic time 
to that mark the "starting time" of the violation, to 
is found to move to higher values as the aging time in- 
creases (Fig. [3Id)). It is interesting to compare to to 
the relaxation times t of the response and the correla- 
tion. We find T fitting the correlation and the response 
with stretched exponential of the form exp[— (i/r)^] and 
(1 — a exp[— (t/r)'^]), respectively. The response ages 
faster than the correlation almost reaching the same re- 
laxation time for the longest i^,- 
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FIG. 3: (a) The interesting region of the FDT plots for different 
aging times: = 90, 125 and 1200 min (cfr. Fig. [2] right panel), 
(b) The inverse slope of the long times points of the FDT plot (see 
Fig . [2J as a function of the aging time, the red line is a guide for the 
eye. This parameter can be interpreted as an effective temperature 
and it is found to reduce to the bath temperature as the aging time 
increases, (c) Waiting time dependence of the Edwards-Anderson 
parameter, the characteristic value of the correlation function at 
which the FDT breaks-down, (d) Evolution with of the char- 
acteristic time at which the FDT is violated (o), together with the 
relaxation time of the correlation function (□) and of the response 
(<0'). The characteristic time for the FDT violation to increases as 
grows followed by the two relaxation times. 



It is important to emphasize that in aU models investi- 
gated so far, for studying the generahzation of the FDR, 
the relaxation time grows roughly as the waiting time: 
T ~ t„ . The aging process that we study experimentally 
here does not obey this simple scaling, the relaxation 
time being several orders of magnitude shorter than the 
typical values of t^. Nevertheless the findings that we 
report in this work indicate that the generalized form of 
theorem applies iitjtw is replaced by t/r in marking the 
transition of the different interesting regions. 

In conclusion, by measuring the autocorrelation func- 
tion of a given variable and the response function of the 
same quantity in an off-equilibrium (aging) colloidal sus- 
pension in the route to the arrested state, we have tested 
the validity of the generalized fluctuation-dissipation re- 
lation. The prediction of the GFDR apply to the present 
experiment, and on the probed time-scale we observe that 
the deviation from the satndard FDT reduces gradually 
as the arrested phase is approached. The characteris- 
tic time at which the violation is seen is always slightly 
above the relaxation time of the measured response and 
correlation function. 
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